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1 Introduction 



Over the past three decades there has been a significant progress in the theory of (1 + 1)- 
dimensional quasihnear systems, 

ul + v]{u)ui = 0, (1) 

which are representable in the Hamiltonian form + P^-'hj = 0. Here h{u) is a Hamiltonian 
density, hj = dy^jh, and P*-' is a Hamiltonian operator of differential-geometric type, 

generated by a metric g^^ (assumed non-degenerate) and its Levi-Civita connection r*j^, via 

6^"' = —g'^^T''gf,. It was demonstrated in [6j that the metric g^^ must necessarily be flat, and 
in the flat coordinates of g'^^ the operator P^^ takes a constant coefficient form P*-' = e*^*-'^- 
In the same coordinates, Hamiltonian systems take a Hessian form ul + e^hijUx = 0. It was 
observed that many particularly important examples arising in applications are diagonalizable, 
that is, reducible to the Riemann invariant form R\ + t;*(R)i?^ = 0. We recall that there exists 
a simple tensor criterion of the diagonalizability for an arbitrary hyperbolic system ([T|) . Let us 
first calculate the Nijenhuis tensor of the matrix u*, 

A/], = v^^duwl - vldy^v] - vi,{dy,vl - dy.v^), (2) 
and introduce the Haantjes tensor 

n% = M;y^vi - M^yy, - Mf.,vy^ + Mf.viv;. (3) 

It was observed in [16] that a (l,l)-tensor Vj with mutually distinct eigenvalues is diagonaliz- 
able if and only if the corresponding Haantjes tensor 7i is identically zero. As demonstrated 
by Tsarev, a combination of the diagonalizability with the Hamiltonian property implies the 
integr ability: all diagonalizable Hamiltonian systems possess an infinity of conservation laws 
and commuting flows, and can be solved by the generalized hodograph transform. We refer to 
[27\ [6] for further discussion and references. 

The aim of our paper is to generalize this approach to (2 + l)-dimensional Hamiltonian 
systems 

ut + A{u)vLx + B{u)viy = 0, (4) 

which are representable in the form ut + P/iu = where h{u) is a Hamiltonian density, and P 
is a two-dimensional Hamiltonian operator of differential-geometric type, 

P^^ = g^^iu)^ + 6'^{u)4 + 9^^{u)^ + 6^(u)n^; 

such operators are generated by a pair of metrics g^^ , g^-' and the corresponding Levi-Civita 
connections r^fc,r^fc via b^j^ = — g^^L-^^, b^^ = — ^**r^^. The theory of multi-dimensional Poisson 
brackets was constructed in [HI \TU\ . The main difference from the one-dimensional situation 
is that, although both metrics g^^ and g'^^ must necessarily be flat, they can no longer be reduced 
to a constant coefficient form simultaneously: there exist obstruction tensors. The obstruction 
tensors are known to vanish if either one of the metrics is positive definite, or a pair of metrics 
is non-singular in the sense of [20], that is, the mutual eigenvalues of g^^ and g^^ are distinct. 
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In both cases, the operator P^^ can be transformed to a constant coefficient form. In the two- 
component situation any non-singular Hamiltonian operator can be cast into a canonical form 



P 



d/dx 
d/dy 



by an appropriate linear change of the independent variables x,y. The corresponding Hamilto- 
nian systems take the form 

uj + (hi)^ = 0, + {h2)y = 0. (5) 

The 'simplest' non-trivial integrable Hamiltonian density is h{u^,u'^) = v}u^ — \{u^Y point 
out that, up to certain natural equivalence, there exist no other integrable densities which are 
polynomial in n^,n^). The corresponding equations ([5]) take the form 

u\ — v}u\ + M^. = 0, uf + Uy = 0, 

see Sect. 4.1. This system appears in the context of the genus zero universal Whitham hierarchy, 
[T71[TH]. Setting = —(pxt, ^ = ^xy one obtains a second order PDE 

1 2 

which is one of the Hirota equations of the dispersionless Toda hierarchy [9j. The same equation 
appeared in [22] in the classification of integrable Egorov's hydrodynamic chains. Other examples 
of integrable Hamiltonian densities expressible in elementary functions include 

etc. The problem of classification of integrable two-component Hamiltonian systems ([5]) was 
first addressed in [10] based on the method of hydrodynamic reductions. We recall that a 
multi-dimensional quasilinear system (jlj) is said to be integrable if it possesses an infinity of re- 
component hydrodynamic reductions parametrized by n arbitrary functions of a single variable 
(see Sect. 2 for more details). It was demonstrated in [TU] that this requirement imposes strong 
restrictions on the corresponding Hamiltonian density h{u^,u^). In Sect. 4 we provide a complete 
list of integrable Hamiltonian densities (Theorem 1), as well as the associated dispersionless Lax 
pairs (Sect. 4.1). The 'generic' density is expressed in terms of the Weierstrass elliptic functions. 
In the three-component situation we consider Hamiltonian operators of the form 



P = 1 — + A2 — , (6) 




here A' are constant and pairwise distinct; the corresponding Hamiltonian systems are 

ui + {hi)x + X\h,)y = 0. (7) 

There is a new phenomenon arising in the multi-component case: it was observed in [12] that the 
necessary condition for integrability of an n-component quasilinear system Q is the vanishing 
of the Haantjes tensor for an arbitrary matrix of the form 

{aA + f3B + -fInyHdA + + 7/„). 
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In fact, it is sufficient to require the vanishing of the Haantjes tensor for a two-parameter fam- 
ily {kA + In)~^{lB + In). We point out that in the two-component case the Haantjes tensor 
vanishes automatically. On the contrary, in the multi-component situation the vanishing of the 
Haantjes tensor is a very strong restriction. Systems with this property will be called 'diag- 
onalizable' (we would like to stress that matrices A and B do not commute in general, and 
cannot be diagonalized simultaneously). In Sect. 5 we obtain a complete list of diagonalizable 
three-component Hamiltonian systems ([7]) (Theorem 3). It turns out that in this case the diag- 
onalizability conditions are very restrictive, and imply the integrability. For technical reasons, 
the classification results take much simpler form when expressed in terms of the Legendre trans- 
form H of the Hamiltonian density h, rather then the Hamiltonian density h itself (recall that 
H = '^u^hi — h, Hi = n*, Uj = hi; we use variables Ui with lower indices for the arguments 
of H). We demonstrate that the Legendre transform H of the 'generic' integrable Hamiltonian 
density h is given by the formula 

A* — A-' 

H = — 7^-iy-V{aiUi, ajUj) 

where 

V{x, y) = Z{x — y) + eZ{x — ey) + e^Z{x — e^y); 

here Oj are arbitrary constants, e = e^'^*/'^, and Z" = (" where (" is the Weierstrass zeta-function: 
C' = — p, (p')^ = ^P^ ~ 53- Notice that we are dealing with an incomplete elliptic curve, g2 = 0, 
and that the expression for V is real. The above formula for H has a natural multi-component 
extension, which is also integrable. This formula possesses a number of remarkable degenerations 
which are listed in Theorems 1 and 3. In particular, one has 

H = — ^ ^ (Q^^^ — CLjUj) ln(ajtij — a-jUj). 

We prove that all examples appearing in the classification possess dispersionless Lax pairs and 
an infinity of hydrodynamic reductions (Theorems 4 and 5 in Sect. 5.1 and 5.2). It is important 
to stress that, in 1 + 1 dimensions, integrable Hamiltonians are parametrized by "^-"^ arbitrary 
functions of two variables. On the contrary, in 2 + 1 dimensions, the moduli spaces of integrable 
Hamiltonians are finite-dimensional. Furthermore, the results Sect. 6 (Theorems 6 and 7) 
make it tempting to conjecture that there exits no non-trivial integrable Hamiltonian systems 
of hydrodynamic type in 3 + 1 dimensions. 

The analysis of the integrability conditions is considerably simplified after a transformation of 
a given Hamiltonian system into the so-called Godunov, or symmetric, form. This construction 
is briefly reviewed in Sect. 3. 

The necessary information on hydrodynamic reductions and dispersionless Lax pairs is sum- 
marized in Sect. 2. 

2 Hydrodynamic reductions and dispersionless Lax pairs 

Applied to a (2 + l)-dimensional system ([3]), the method of hydrodynamic reductions consists 
of seeking multi-phase solutions in the form 

u(x, y, t) = u{R\x, y, t), i?"(x, y, t)) 
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where the 'phases' R^{x,y,t) are required to satisfy a pair of (1 + l)-dimensional systems of 
hydrodynamic type, 

Rl = J^*(R) Rl, Rl = At*(R) R'y 

Solutions of this form, known as 'non-hnear interactions of n planar simple waves' |25l [U I24j . 
have been extensively discussed in gas dynamics; later, they reappeared in the context of the 
dispersionless KP hierarchy, see \14:\ [T5] and references therein. Technically, one 'decouples' 
a (2 + l)-dimensional system (j3]) into a pair of commuting n-component (1 + l)-dimensional 
systems. Substituting the ansatz u(i?^, ...,R^) into (jH) one obtains 

{u'ln + fi'A + B) diu = 0, i = l,...,n, (8) 

di = d/dR^, implying that both characteristic speeds z^* and /i* satisfy the dispersion relation 

det(i//„ + fiA + B) = 0, (9) 

which defines an algebraic curve of degree n on the (v, /i)-plane. Moreover, and have to 
satisfy the commutativity conditions 

= Jil^, (10) 

— — /i* 

i ^ j, see |27) . In was observed in [TU] that the requirement of the existence of 'sufficiently 
many' hydrodynamic reductions imposes strong restrictions on the system and provides an 
efficient classification criterion. To be precise, we will call a system @ integrable if, for any n, 
it possesses infinitely many n-component hydrodynamic reductions parametrized by n arbitrary 
functions of a single variable. Thus, integrable systems are required to possess an infinity of 
n-phase solutions which can be viewed as natural dispersionless analogs of algebro-geometric 
solutions of soliton equations. 

We recall that a system (jl]) is said to possess a dispersionless Lax pair 

i^t = f {u,i)y) , V'x = 5 (u, -02/) , (11) 

if it can be recovered from the consistency condition ipxt = ''Ptx (we point out that the dependence 
of / and g on ipy is generally non- linear). Lax pairs of this type first appeared in the construction 
of the universal Whitham hierarchy, see [17] and references therein. It was observed in [28] that 
such non-linear Lax pairs arise from the usual 'solitonic' Lax pairs in the dispersionless limit, and 
the cases of polynomial/rational dependence of / and g on ■0^ were investigated. In particular, 
a Hamiltonian formulation of such systems was uncovered, requiring a non-local Hamiltonian 
density. It was demonstrated in [101 113j that, for a number of particularly interesting classes of 
systems, the existence of a dispersionless Lax pair is equivalent to the existence of hydrodynamic 
reductions and, thus, to the integrability. 

Setting Ipy = p and calculating the consistency condition ipxt = iptx by virtue of ([!]), one 
arrives at the following relations for f(u,p) and g{u,p): 

grad/ + grad^ ^ = 0, grad g [fpin +gpA + B] = 0; (12) 

here grad is the gradient with respect to u. In particular, this shows that fp and gp satisfy the 
dispersion relation Q , and the vector grad g belongs to the left characteristic cone of the system 
Thus, as p varies, the equations u = fp, fJ, = gp parametrize the dispersion curve ([9]), while 
grad g parametrizes the left characteristic cone. 

Throughout this paper we assume that the dispersion relation ([9]) defines an irreducible 
algebraic curve. This condition is satisfied for most examples discussed in the literature so far. 
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3 Transformation of a Hamiltonian system into Godunov's form 



Recall that a system of hydro dynamic type ^ is said to be symmetrizable, or reducible to 
Godunov's form [8], if it possesses a conservative representation of the form 

{du-p)t + {du^q)x + {du-r)y = 0; 

here the potentials p, q and r are certain functions of u. Any such system possesses an extra 
conservation law L{p)t + L{q)x + L{r)y = where L denotes Legendre's transform. Equations 
in Godunov's form play important role in the general theory of multi-dimensional hyperbolic 
conservation laws [5j. 

Given a Hamiltonian system ([7]) we perform the Legendre transform, H = L[h) = u^hi — 
h, Hi = u^, Ui = hi, to obtain a system in Godunov's form, 

iHi)t + iui)^ + X\ui)y = 0, 

which corresponds to the choice p = H, q = Yl"^! /'^^ ^ — S-^*'"?/^- We assume that the 
Legendre transform is well-defined, that is, all partial derivatives hi are functionally independent. 
This condition is equivalent to the requirement that the Hessian matrix of h is non-degenerate, 
which is automatically satisfied under the assumption of the irreducibility of the dispersion 
relation. It turns out that the integrability conditions take much simpler form when represented 
in terms of the Legendre transform H = L{h), rather then the Hamiltonian density h itself. Thus, 
in what follows we will work with systems represented in Godunov's form (to make the equations 
look formally 'evolutionary' we will relabel the independent variables as x,y,t T, X, Y). This 
results in 

{ui)T + X\ui)x + {Hi)Y = 0; (13) 

Systems of this type can be viewed as describing n linear waves (traveling with constant speeds 
A* in the X, T-plane) which are non-linearly coupled in the y-direction. 

4 Integrable Hamiltonians in 2 + 1 dimensions: two-component 

case 

In this section we classify two-component Hamiltonian systems The corresponding Legendre 
transform is 

VT + {Hy)Y = 0, WX + {H^)y = 0; (14) 

here v = ui,w = U2- We point out that this case was addressed previously in [10], although the 
classification was only sketched. Here we provide a complete list of integrable potentials H{v, w), 
and calculate the corresponding dispersionless Lax pairs. For systems (jl4p the integrability 
conditions constitute an over-determined system of fourth order PDEs for the potential H{v, w): 



(15) 
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The system ([13]) is in involution, and its solution space is 10-dimensional [lOj. We point out 
that the transformations 

V ^ av + b, w ^ cw + d, H ^ aH + + 711;^ + fiv + uw + 5 

generate a 10-dimensional group of Lie-point symmetries of the system (jlSp . These transforma- 
tions correspond to obvious linear changes of the independent variables X, Y, T in the equations 
(I14p . One can show that the action of the symmetry group on the moduli space of solutions of 
the system (jlSp possesses an open orbit. The classification of integrable potentials H{v, w) will 
be performed up to this equivalence. Moreover, we will not be interested in the potentials which 
are either quadratic in w and generate linear systems (|14p . or separable potentials of the form 
f{v) + g{w) giving rise to reducible systems. Our main result is the following complete list of 
integrable potentials: 



Theorem 1 The 'generic ' solution of the system il5\) is given by the formula 

H{v, w) = Z{v + w) + eZ{v + ew) + e^Z{v + e'^w); 

here e = e^'^*/^ and Z"(s) = C{s) where C is the Weierstrass zeta- function: 
—53. Degenerations of this solution correspond to 

H{v,w) = ^v^C{w), 

H{v, w) = {v + w) ln(v + w), 
as well as the following polynomial potentials: 

H{v, w) 



2 2 
V w , 



H{v, w) 



a 



vw H — w 
5 



a = const, 



and 



H{v, w) 



vw H — w . 
6 



(16) 

-P, ip'? = 

(17) 
(18) 

(19) 
(20) 

(21) 



Remark. The 'elliptic' examples (1160 and ()17p possess a specialization g^ 
({"w) — > l/ui, a{w) — > w, etc. This results in the potentials 



0: p{w) 



H{v, w) = {v + w) log{v + w) + e{v + ew) log{v + ew) + e^(f + e'^w) log{v + e^w) 



(22) 



and 



H{v,w) = -, 



respectively. Dispersionless Lax pairs for the equations ([T 
(fTH]) - (pT|) are calculated in Sect. 4.1. 

Proof of Theorem 1: 



corresponding to the potentials 



vw 



The system ([TC]) can be solved as follows. The first two equations imply that H^^/H^ 
const. Similarly, the last two equations imply H^ww/H^w — const. Setting H^^ = e one can 
parametrise the third order derivatives of H in the form 



1 



-me 



H„ 



1 



-ne 



(23) 
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here m, n are arbitrary constants. The compatibihty conditions for these equations, plus the 
equation ([ISjs, result in the following overdetermined system for e: 

, , TY}jTh 2 /o/i\ 

The general solution of the first (Liouville) equation has the form 

g2 ^ 4 p'{v)q'{w) 

mn {p{v) + q{w)y ' 

one has to consider separately the case e = const (up to the equivalence transformations, this 
results in the potential ([2T]) ). as well as the case when e depends on one variable only, say, on w 
(this leads to the potential ()20|) ). Let us assume that both constants m and n are nonzero (the 
cases when either of them vanishes will be discussed later). By scaling v and w one can assume 
m = n = 1. Setting 

{p'f = P\p), {q'f = Q\q), (26) 

(here P{p) and Q{q) are functions to be determined), one obtains from the last two equations 
()24p the following functional-differential equations for P and Q: 



P"{p + qf-4P'{p + q) + 6P = 2Q'{p + q)-6Q, Q"{p + qf -4Q'{p + q) + 6Q = 2P'{p + q)-6P; 

these equations imply that both P and Q are cubic polynomials in p and q, 

P = ap^ + bp^ + cp + d, Q = aq^ - hq^ + cq - d, 

where a, b, c, d are arbitrary constants. Notice that the right hand side of possesses the 
following 5L(2, i?)-invariance, 

ap + [5 ap — P 

7P + jp — 

which can be used to bring the polynomials P{p) and Q{q) to canonical forms. There are three 
cases to consider. 

Three distinct roots: in this case one can reduce both P{p) and Q{q) to quadratics, so that 
the ODEs ()26p assume the form 

ipr = Y^' + 93? and iqr = Y(^'+93f, 

respectively. Thus, p = p'{v), q = p'{w) where p is the Weierstrass p-function: (p')^ = 4p^ — gs 
(we point out that the value of is not really essential, and can be normalized to ±1). Setting 

rr _ _ 12p(7;)p(u;) 



p'{v) + p'{w) 

and integrating (I23p with respect to v and w we obtain 



H^^--bQ{w) — — J-—, — — — — , Hu,y, - -b(;[v) —— (27) 

p'{v) + p'{w) p'{v) + p'{w) p'[v) + p'[w) 
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here the zeta- function is defined as C = —p. Since the p- function on the elhptic curve = 
— satisfies the automorphic property p(ez) = ep(z), e'^ = 1, one can rewrite (f27|) in the 
following equivalent form: 



Hvv = -2\Q(v + w) + eC{v + ew) + eX{v + ew 
Hvw = -2{c,{v + w)+ e^C{v + ew) + eC{v + e^w 
Hu,u, = -2U{v + w) + C{v + ew) + C{v + e^w) 



Up to a constant multiple, these formulae give rise to (jl6|) . 

Double root: in this case both P{p) and Q{q) can be reduced to p and q, so that the ODEs 
(|26p take the form {p')'^ = 27p^ and {q')^ = 27g^, respectively. This leads io p = , q = w'^, 
and a straightforward integration of ()23p gives 



6w 6vw 6v 

_|_ yj^ yO _|_ y^O _|_ yjO 

notice that these formulae can be obtained as a degeneration of (|27p corresponding to = 0. 
Up to a constant multiple, this leads to the potential ()22p . 

Triple root: in this case both P{p) and (5(g) can be reduced to constants, so that the ODEs 
(f26l) take the form {p')^ = 1 and {q')^ = 1, respectively. This leads to e = 2/(t; + w), which, up 
to a constant multiple, results in the potential (jlSp . 

If m = 0, n 7^ (without any loss of generality we will again set n = 1), equations (|24p can 
be solved in the form e = 6vp{w) where p is the Weierstrass p-function: (p')^ = 4p^ — 33. The 
corresponding potential H is given hy H = —3v'^C{w). Up to a multiple, this is the case p!7|) . 

In the simplest case m = n = equations ()24|) imply 

e = {av + (3){'yw + 5), 
and the elementary integration of equations p3p results in 

H{v, w) = i^oiv'^ + Pv){^^w'^ + 6w); 

here a, f3, 7, (5 are arbitrary constants. Using the equivalence transformations one can reduce H 
to either H = v^w"^ (both a and 7 are nonzero) or H = vw"^ (a = 0). These are the polynomial 
cases (JTH]) and a subcase of (pUp . respectively. This finishes the proof of Theorem 1. 



4.1 Dispersionless Lax pairs 

In this section we calculate dispersionless Lax pairs for systems (|14p corresponding to the po- 
tentials (fT6]) - ([2T]) of Theorem 1. We point out that, in spite of the deceptive simplicity of some 
of these potentials, the corresponding Lax pairs are quite non-trivial. 

Potential (12 ip : The corresponding system ()14p takes the form 

vt + wy = 0, wx + wwy + = 0; (28) 

it arises in the genus zero case of the universal Whitham hierarchy [17 t ll8|. This system possesses 
the Lax pair 

i;T = ^H'^Y + w/2), i)x=^l + v/2. 



9 



A simple calculation shows that the Legendre transform of the potential H{v,w) = vw + ^w'^, 
defined by the formulae 

= Hy, = H^, h{v} ,u^) = vH^ + wH^ — h, 

is also polynomial: 

6 

We point out that all other examples of integrable potentials H{v, w) produce non-polynomial 
Hamiltonian densities h{u^,u'^). 

Potential (I20p : The corresponding system ()14p takes the form 

VT + {w'^)y = 0, wx + 2{vw)y + a{w^)Y = 0. (29) 
For Q = it possesses the Lax pair 

V'T = -7772-, tpx = tpY - '^Vlpy . 

Setting V = uy, w"^ = —uj- one can rewrite (|29p (when a = 0) as a single second order PDE 

UXT + 2uyuty + ^utuyy = 0. 

Up to a rescaling X — > —2X this equation is a particular case of the generalized dispersionless 
Harry Dym equation [H |23]. For a 7^ the Lax pair modifies to 

'^T = f (^^) ' V'x = V'y - ^VIpY, 

where the function f{s) satisfies the equation f'{s) = —s/{as^ + 1) (for a = one recovers the 
previous formula). The first equation of this Lax pair appeared in [23| as a generating function 
of conservation laws for the Kupershmidt hydrodynamic chain. Without any loss of generality 
one can set a = — 1, which gives 

/(s) = \ (ln(s - 1) + ln(s - e) + e ln{s - e^)) , ^ 1. 

Potential (I19|) : The corresponding system (|14p takes the form 

VT + 2{vw^)y = 0, WX + 2{v^w)y = 0. (30) 

It possesses the Lax pair 

where the dependence of a and b on tpy = ^ \s governed by the ODEs 

a' = -4? - 2, 6' = 4- + 2. 
a 

To solve these equations we proceed as follows. Expressing b from the first equation, b = 
—A:a/{a' + 2), and substituting into the second one arrives at a second order ODE 2aa" — 
3(a')^ + 12 = 0. It can be integrated once, (a')^ = Aca? + 4, where c is a constant of integration. 
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Without any loss of generality we will set c = 1. Thus, a is the Weierstrass p-function: a = 
p(^, 0, —4) = p(^). The corresponding b is given by 6 = — 4p/ (p'+2). Notice that this expression 
for b equals p(^ + c) where c is the zero of p-function such that p(c) = 0, p'(c) = 2 (use the 
addition theorem to calculate p(^ + c)). Ultimately, we obtain the Lax pair 

Setting V = v'^, W = up' one can rewrite (I30p in the form where the non-linearity is quadratic: 



Vt + 2WVy + 4yVFy = 0, + 2VWy + 4M/yy = 0. 



'^T = -f{w,tpY), tpx = --V b{ip- 



Y 



Potential ( II Sh : The corresponding system (I14p takes the form 

VY + WY „ VY + WY „ 

H \ — = 0, wx^ ■ — = 0. 

V +W V + w 

It possesses the Lax pair 

V'T = — ln(t(; -|- V'y), 4'x = — V'y)- 

This system also arises in the genus zero case of the universal Whitham hierarchy |17[ dB] ; its 
dispersionful analogue was constructed in [26j . 

Potential (II 7p : The corresponding system (|14p takes the form 

Vt + C,{w)vy — vp{w)wY = 0, wx — p{w)vvy — -v'^p'{w)wy = 0. 
One can show that it possesses the Lax pair 

1 

— '( 
2 

where, setting ^pY = the function /(tf,^) has to satisfy the equations 

_ 2b{0p{w) f . , VH 
b'iO + p'H' ^^""'^ b'{0 + p'{w)- 

We point out that the consistency condition f^^ = f^^ implies a second order ODE 2bb" — 
3{b'f - 353 = which, upon integration, gives 

(6'(0)' = 463(0 -53, 

(the constant of integration is not essential). Thus, one can set b = p{^) so that the equations 
for / take the form 

Jp{OpH_ 2p\w) 

P'iO + p'H' ^ P'(0 + p'H' 

compare with ([27|) ! Thus, 

f{w,0 = - \na{^ + w) + - lna{^ + ew) + jlna{^ + e^w), 
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where a is the Weierstrass sigma-function: a' /a = C- Ultimately, the Lax pair takes the form 

ipT = -lnfT(V'y + w) + - In a{tpY + ew) + y lncj(V'j/ + e^w), ipx = --v^p{ipY)- 
Potential (I16p : the equations corresponding to H/?, take the form 

One can show that the corresponding Lax pair is given by the equations 

i^T = f{w,'lpY), ^X=9{v,tl^Y) 

where, setting ■^y = ^, the first order partial derivatives of / and g are given by 

2p{0p{w) 2p\w) 



P'(0 + P'H' ' ' ' p'{i) + p'{w) 
and 

respectively. Explicitly, one has 

f{w,0 = -3ln(T(C + u;) - - lna{^, + ew) - y lncr(^ + e^w), 

9{v,0 = 3 lncr(^ - v) + -lncr(C - ev) + —lna{^ - e^v). 

Notice that the expression for f{w, ^) coincides with the one from the previous case. This means 
that the corresponding Hamiltonian systems commute with each other — the fact which is, in 
a sense, unexpected. 

Potential ( 122^ : this is the 53 = degeneration of the potential (I16p . The system corresponding 
to H/3 takes the form 

w'^ vw vw 

'^T + — 3 — — '^'^y = 0' — — '^'"y + — ^'^y = 

_|_ y^J -yd _j_ yj^ y6 _j_ yj6 _j_ y^J 

it possesses the Lax pair 

V'T = /{w/tPy), ipx = giv/ipY) 

where the dependence of / and g on their arguments is specified by f'{s) = s/{s'^ — 1), g'{s) = 
s/{s^ + 1). Explicitly, one has 

fis) = ^ {Hs - 1) + ln(s - e) + eln(s - e^)) , 
g{s) = (ln(s + 1) + e2ln(s + e) +eln(s + e2)) . 
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5 Integrable Hamiltonians in 2 + 1 dimensions: three-component 
case 



In this section we classify three-component integrable equations of the form (1131) . 

= 0, (31) 








^ 
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H23 


Hss ) 
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assuming that the constants A* are pairwise distinct. As mentioned in the introduction, the 
integrability of the system (I3ip implies the vanishing of the Haantjes tensor for any matrix of 
the two-parameter family {kA-\- Is)~^{lB + Is). Here A = diagi^X) and B = {H-ij). To formulate 
the integrability conditions in a compact form we introduce the following notation: 

Ri = ^^[X^-X'), R2 = ^^{X'-X'), Rs = ^^{X'-X'); 

J^23 ^13 ^12 

we will see below that all mixed partial derivatives Hij must be non-zero, otherwise the system 
is either linear, or reducible. Moreover, we will need the quantities 

/ = - 4(A2 - X^){X^ - X^)H^2 - 4(A^ - X^){X^ - X^)hI^ - 4(A^ - X^)[X^ - X^)Hf^ 

and 

J = (A^ — A^)iJi2-f^i3 + i^^ ~ ^^)H2sHi2 + (A^ — X^)Hf^H2s — Hi2H2sHisA 

where 

A = (A2 - X^)Hii + (A^ - X^)H22 + (Ai - X^)Hss. 
Our first result is the following 



Theorem 2 The system ^31]) with an irreducible dispersion curve is diagonalizable if and only 
if the potential H satisfies the relations 

J = 0, Hi2s = 0, 
^ ((A3 - X^)Hii +R2 + Rs) = 0, 



(32) 

^((Al-A3)F22+iil + i?3) =0, 



^((A2-Al)F33+iil + i?2) =0. 

Notice that, in contrast to the two-component situation (llSp . these relations are third order in 
the derivatives of H. We will demonstrate below that the necessary conditions (|32|) are, in fact, 
sufficient for the integrability, and imply the existence of dispersionless Lax pairs and an infinity 
of hydro dynamic reductions. 

Remark. The condition J = 0, which is equivalent to Ri + R2 + Rs = A, has a simple 
geometric interpretation as the condition of reducibility of the left characteristic cone of the 
system (I3ip (see Sect. 2 for definitions). Indeed, the left characteristic cone consists of all 
vectors g = (91,02,93) which satisfy the relation 

g(z./3 + fiA + B) = 0. (33) 
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Excluding v and ^, one obtains a single algebraic relation among gi,g2,g3, 
{Hi3{gifg2 + H2zgi{g2? + i^335i5253) (A^ - \^) + 

{H2i{g2fg3 + ^1352(53)' + i^ii5i5253) (A' - X^)+ (34) 
(^23(53)^51 + Hugaigif + ^^22515253) (A^ - A^) = 0, 

which is the equation of the left characteristic cone. The condition J = is equivalent to its 
degeneration into a line and a conic: 

[HuHisgi + Hi2H23g2 + -f^i3-f^235'3] , , 

[Hi3H23{X' - X^)gig2 + ^i2i^23(A=' - X^)gm + H^HisiX^ - A=')5253] =0. 

We point out that, by virtue of (I33p . the left characteristic cone and the dispersion curve are 
birationally equivalent. This implies that the dispersion curve is necessarily rational, although 
not reducible (the linear factor of the left characteristic cone corresponds to a singular point on 
the dispersion curve — see Sect. 5.2 for explicit formulae). 

Proof of Theorem [2} 

To simplify the calculation of the Haantjes tensor we multiply the matrix {kA+l3)^^{lB-\-I^) 
by {kX^ + l){kX'^ + l){kX^ + 1). This results in the matrix A{IB + 13) where A = diag[{kX^ + 
1)(/cA3 + 1), {kX^ + l){kX^ + l), {kX^ + l){kX^ + 1)]. Since the multiplication by a scalar does not 
effect the vanishing of the Haantjes tensor, we will work with the matrix A[IB + L^) which has 
an advantage of being polynomial in k and I. Using computer algebra we calculate components 
of the Haantjes tensor Ti (which are certain polynomials in k and /) and set them equal to zero. 
First of all, one can verify that all components of the form 7i\j vanish identically, so that the 
only nonzero components are T~L^jk, i ^ j ^ k. In the following we will focus on the analysis of 
the component ^^2- it turns out the vanishing of alone implies the vanishing of the full 
Haantjes tensor. Let us compute coefficients at different powers of the parameter / and set them 
equal to zero. At the order , all terms in TL12 vanish identically since A is a constant diagonal 
matrix. The coefficient at is a polynomial in /c, however, setting its coefficients equal to zero 
we obtain only one independent relation: 

-f^l23 = 0. 

Similarly, two extra relations come from the analysis of /^-terms, three relations from /^-terms, 
and four relations from /^-terms. Ultimately, we end up with a set of 9 linear homogeneous 
equations for the 9 third order derivatives Hm, Hnj. 

From these 9 relations it readily follows that if one of the mixed derivatives equals zero, say, 
H12 = 0, then either H13H23 = or Hij^ = for all i,j,k. In the first case the system ()3ip 
decouples into a pair of independent 1x1 and 2x2 subsystems. The second case corresponds 
to linear systems with constant coefficients. Therefore, from now on we assume Hij ^ for any 

The set of 9 relations so obtained is rather complicated, and the calculation of the corre- 
sponding 9x9 determinant is computationally intense. A simpler equivalent set of relations 
can be derived as follows: first, divide TI12 by {X^k + 1)(A^A: + l)^(A^/c + 1)^ (which is a com- 
mon multiple), then equate to zero the coefficient of l'^ at k = — 1/A^,— 1/A^ (the coefficient 
at k = — 1/A^ appears to be a linear combination of the previous two), the coefficient of at 
k = -l/A^, -1/A^ -l/A^ and the coefficient of l'^ at k = 0, -1/A\ -l/A^, -l/A^. As a result. 
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we arrive at a simpler set of 9 linearly independent relations that are nothing but linear combi- 
nations of the previous ones. If the determinant of this system is non-zero, then all remaining 
derivatives Hm and Haj vanish identically. This is the case of linear systems. Thus, to obtain 
non- linear examples, one has to require the vanishing of the determinant. It is straightforward 
to verify that this determinant factorizes as follows: 

Thus, there are two cases to consider. If 

/2 - 64(Ai - A2)(a2 - a3)(A3 - A^) J = 0, (36) 

then the dispersion relation of the system (j3ip is reducible. To show this we introduce the 
quantities 

01 = AFi2-2Fi3/723(A^-A2), 

02 = AF23 - 2Hi2Hi^{\^ - A^), 

03 = A2-4i/23(Ai-A2)(A2-A3), 

04 = hU)? - A2) + hU\^ - A^), 

which can be verified to satisfy the quadratic identity 

{x^ + -\^)^l + v^^^li = {). (37) 

In terms of these quantities, the equation (|36p can be rewritten as follows: 

(O3 - 4(Ai - \^)aif + 16(Ai - \^){\^ - = 0, (38) 

or, equivalently, 

(O3 + 4(a1 - \^)^if + 16(A^ - \^){\^ - >?)al = 0; (39) 

one has to use the identity (f37|) to verify the equivalence of ([381) and ([391) . Let us assume that 
A^ < A^ < A^. Since we are interested in real- valued solutions, the equation (|38|) implies 

= 0, 1^3 = 4(A^ - A^)174; (40) 
(one should use ([39]) if A^ < A^ < A'^ ) . In this case the identity ([37|) takes the form 

(A^ - \^)^^{ + 4(A^ - = 0, 

so that 

f]^ = 0, il4 = 0. 
These conditions lead to potentials of the form 

H = U2{'yui + dus) + f{jui + 6U3); 

here the constants 7 and 6 satisfy the relation (A^ — A"'^)^^ -|- (A^ — A'^)7^ = 0, and / is an 
arbitrary function of the indicated argument. This ansatz, however, implies the reducibility of 
the dispersion relation as discussed in [12]. Thus, we are left with the second branch J = 0, 



15 



in which case the rank of the system drops to 5, and we end up with the equations (|32|). This 
finishes the proof of Theorem 2. 

The main result of this Section is a complete list of integrable potentials i7(Mi, M2, M3) which 
come from a detailed analysis of the equations (I32p . The classification will be performed up to 
the following equivalence transformations, which constitute a group of point symmetries of the 
relations (f32]) . 

Equivalence transformations: 

transformations of the variables Ui'. Ui ^ aui + hi] 
transformations of the potential H: 

the latter corresponding to y — > oY + [3T + jX in the equations ()3ip . Moreover, relations 
(j32p are invariant under arbitrary permutations of indices. Finally, we will not be interested in 
the potentials which are either quadratic in Ui and generate linear systems (I3ip . or separable 
potentials, e.g., H = f{ui) + g(u2,U3), giving rise to reducible systems. 



Theorem 3 The 'generic ' solution of the equations ^3B^ is given by the formula 

A* — A-' 

-H' = -X 2 2 V{aiUi,ajUj) (41) 

where 

V{x, y) = Z{x -y) + eZ{x - ey) + e^Z{x - e^y); (42) 

here e = e^'^*/^ and Z" = ( where ( is the Weierstrass zeta- function: (' = —p, (p')^ = —93- 
Degenerations of this solution correspond to 

A* — A-' 

-ff = - X 22 V{aiUi,ajUj) (43) 
■ / - 3(X,- a^ 

where 

y(x, y) = {x — y) ln(x — y) + e(x — ey) ln(x — ey) + e^{x — e^y) ln(x — e^y), 

and 

_^ yi _ 

H = — X^ — 2 2 ("»^» ~ ^i^i) ~ ^j^j)^ (^4) 

respectively. Further examples include 

H = J^ulC{a2U2) + 2~-'"lC(«3'"3) - 2 2 ^(«2^2,«3^3) (45) 

"2 ^3 "2 "^3 

where V is the same as in ^J^- This potential possesses a degeneration 

H={X'- X^)ujul + (A2 - X^)C{us + c)ul - (A3 - X')C{ns)ul (46) 

here C,' = — p, (p')^ = ^p^ + 4, and c is the zero of p such that p(c) = 0, p'(c) = 2. possesses 
a further quartic degeneration, 

H = {X^ - X'^)ulul + (A^ - >^^)ulul + (A^ - A^)ui?x?. (47) 
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We have also found the following (non- symmetric) examples: 

H = {pui + qu's) In {pui + qu3) - gp(A^ - X'^){X^ - A^)'U,f- 
ig(A3 - a1)(A3 - X^)ul + p{X^ - X^)u2U3 + q{X^ - X^)uiU2, 



(48) 



H = {X^- X^)u2ul + (A^ - X^)u2ul + — (A^ - a3)(a3 - X^)ul + ^, (49) 

iU W3 



and 



(50) 



where 



H = {X^- X^)u2ul + (A2 - A3)u2txi+ 
^(A2 - Ai)(Ai - A3)^f + ^(A2 - A3)(A3 - Ai)^ + u,G(p}^, 

G{x) = (px + q) log (px + g) + e(px + eq) log (px + eq) + e^(px + e^g) log {px + e^q). 

Up to the equivalence transformations, the above examples exhaust the list of integrable potentials. 
We claim that all examples appearing in the classification possess dispersionless Lax pairs and 
an infinity of hydrodynamic reductions (this will be demonstrated in Sect. 5.1-5.2). 

Proof of Theorem [3} 

We can assume that all mixed partial derivatives Hij are non-zero. It follows from (|32p that 

f Hi2His\ 9^ f Hi2H23\ f HisH23 i „ 

0. (51) 



duidu2du3 \ H23 J duidu2du3 \ H13 J duidu2du3 \ H12 
The further analysis depends on the value of the expression 

du2 du3 dui dui du2 du^ ' 



which appears as a denominator when solving the equations (j5ip . 

Case I. The expression (I52p is nonzero. In this case equations ()5ip are equivalent to 



F K F K 

*-^M2,M3 — p '^uz ~r ^U2 F K 

J^U3,ui — Q i" p -f^-na G F 



where F = I/H12, G = I/H23, K = I/H13. Keeping in mind that F3 = Gi = K2 = 0, we can 
rewrite these equations in the form 

^-«-^--»- <^') 



The system (|53P possesses obvious symmetries 

F ^ fl{ui)f2{u2)F, G^ f2{u2)f3{u3)G, K ^ fi{ui)f3{u3)K, 

(54) 

ui^gi{ui), U2^g2{u2), U3^g3{u3); 
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here and gi are six arbitrary functions of the indicated arguments. As a first step, we introduce 
the new variables 

^ K F' F G' G K' 

which are nothing but the invariants of the first 'half of the symmetry group (|54p . In terms of 
p, q, r, the equations ([53]) take the form 



qi = -P2 = pq, r2 = -qs = qr, ps = -ri = pr. (55) 

This system is straightforward to solve: assuming p ^ (the case when p = q = r = will be 
a particular case of the general formula), one has q = —p2/p, f = P^/p-, along with the three 
commuting Monge- Ampere equations for p, 

P23 = 0, (lnp)i2 =P2, (lnp)i3 = -P3. (56) 

The integration of the last two equations implies pi/p = p + 2ip{ui,U3) and pi/p = —p + 
2ip{ui,U2), respectively. Thus, p = ip{ui,U2) — ip{ui,U3), and the substitution back into the 
above equations gives -01(^1) ""2) — '^^(^1)^2) = ^liui^u^) — ip'^{ui,U3). The separation of 
variables provides a pair of Riccati equations, ipi = ip'^ + V{ui) and ifi=if'^ + V{ui). Thus, 
■0 = — [Infji, ip = — [Inwji, where v and v are two arbitrary solutions of the linear ODE 
vii + V{ui)v = 0. Therefore, we can represent ip and 99 in the form 

= -[In{q2{u2)pi{ui) - P2iu2)qi{ui))]i, ip = -[In(g3(n3)pi(ui) - qi{ui)p2{u2))]i, 

where pi{ui) and qi{u\) form a basis of solutions of the linear ODE. Introducing Wi{ui) = 
qi{ui)/pi{ui), one obtains the final formula 

W[{W3-W2) 

P = Tp — if — 



{W2 - Wi){w3 - Wi) ' 

leading to 

W2{wi — W3) w'^{w'^ — W^) 



{W2 - Wi){w2 - {W3 - Wi){w3 - W2)' 

Here Wi{ui) can be viewed as three arbitrary functions of one argument. The corresponding 
F,G, H are given by 

F = SiS2iwi - W2), G = 8283(102 - W3), K = 81S3{W3 - Wi), 

where Si{ui) are three extra arbitrary functions. This implies the ansatz 

_ P{UI)Q{U2) ^ _ Q{U2)R{U3) ^ _ P{UI)R{U3) ^ 

fiui)-g{u2)' g{u2)-h{u3y h{u3)-f{ui)' 

(with the obvious identification wi{ui) f{ui), si(ni) 1/P{ui), etc). We have to consider 
different cases depending on how many functions among /, g, h are constant. 
Subcase 1: f = g' = h' = 0. Without any loss of generality one can assume 

H12 = P{UI)Q{U2), H23 = Q{U2)R{U3), Hi3 = Piui)R{u3). 



Substituting this ansatz into (|32p one can show that the functions P, Q, R must necessarily be 
linear. Up to the equivalence transformations, this leads to a unique quartic potential (I47p : 

H = {X^ - X'^)ulul + (A^ - X^)ulul + (A^ - X^)ulul 
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Subcase 2: f = g' = 0. Without any loss of generality one can assume the following ansatz: 

Hl2=P{ui)Qiu2), H23 = Qiu2)hi{u3), i/13 = )/i2 (^^3) . (58) 

The substitution into (j32p implies that P and Q must necessarily be linear. Up to the equivalence 
transformations, this results in the potential 

where the functions a and b satisfy the ODEs 

a" = A-2, b" = A^-2, a'b' = 2{a + b). 
b' a' 

The special case a = b = u\ brings us back to the quartic potential from the previous subcase. 
The generic solution of these ODEs takes the form a{u3) = —Ci^s), b{u3) = C{u3 + c) where 
C iz the Weierstrass (^-function, ^' = — p, (p')'^ — ^P^ + 4, and c is the zero of p such that 
p(c) = 0, p'(c) = 2. This is the case (li^ . 

Subcase 3: /' = 0. The analysis of this case leads to the ansatz 

H = {X^- X^)ula{u2) + {X^ - X^)ulb{u3) + h{u2, ua) 
where ^ ^ 

(here 02,03 are arbitrary constants), and the second order derivatives of h{u2,U3) are given by 

A^-A^ p(o2t^2)p(a3U3) 
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0203 p'(a2n2) - p'(a3n3) ' 



rr , A^-A3 fl^^ ^ p2(03U3) 

H22 = 4 2 — 0^(03^3^ 



03 \2 p' {a2U2) - p' [asus] 

P''[a2U2) 



H33 = 4 2 — oC(a2tt2) . 

0.2 V2 P (03^3) - P (a2^t2 

This is the case (j45|) . 

Generic subcase: /'(x) ^'(x) /i'(2;) 7^ 0. From ()57p and ()32|) we find all third order derivatives of 
H. The compatibility conditions diHjji = djHuj give rise to six functional-differential equations 
for the functions f,g,h,P,Q,R. It follows from ()32p that 



+ (Ai - A3)i722 + (A2 - Al)i733j = 0, 
du, (R2 + (A2 - X^)H33 + (A3 - A2)ii'n) = 0, (59) 

dm {R3 + {X^ - >?)Hii + (Ai - X^)H22) = 0. 

These give us three more equations for /, h, P, Q, R, so that we have nine equations altogether. 
Substituting the values of the third order derivatives of H into the first equation ()59p . taking 
the numerator and dividing by the common factor P{ui)^Q{u2)'^R{u3), we get a fourth degree 
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polynomial in f,g,h,P,Q,R, and first order derivatives thereof. Applying to this polynomial 
the differential operator 

/'(ni) g'{u2) f'{ui)g'{u2) h'{u^) 
we arrive at a separation of variables, 

(A^ - \^){P"'{u^)nu^) - P"{u,)f"im)) (A3 - X^)iQ"'iu2)g'iu2) - Q" {u2)g" {u^)) 

Integrating twice, we obtain 

(A2 - a3)P' = cf + aif + h, (A^ - X^)Q' = cg^ + 035 + 63. 

Analogously, 

(A^ -A^)/?' = 0/1^ + 035 + 63- 

Using these relations we eliminate all derivatives of P, Q and R from our nine equations. As a 
result, we obtain a linear system of nine equations for the three unknowns P, Q, R. This system 
is consistent (that is, the rank of the extended matrix is < 3) if and only if ai = a2 = 0.3 = a, 
bi = b2 = = b, and 

4(c/i2 + 6/1 + a)h'^f" - 4(c/2 + bf + a)f'^h"+ 
if - h) f 2c(/2 + fh + h^) + 36(/ + h) + 6a) f'h" = 0, 



4(c/2 + bf + a)f'^g" - A{cg^ + bg + a)g'^f"+ 
(g-f) (2c{g^ + gf + f) + 3b{g + /) + 6a) 5"/" = 0, 

4(cff2 + bg + a)g'^h" - 4{ch^ + bh + a)h'^g"+ 
[h - g) f 2c(/i2 + hg + g"^) + 3b{h + g) + 6a) h"g" = 0. 



(60) 



Hence, we have either c = b = a = or f" = g" = h" = 0, otherwise f"g"h" 7^ 0. If 
c = 6 = a = then the linear system for P,Q,R becomes homogeneous. Its rank equals two if 
and only if /" = g" = h" = 0. In this case 

Pf'{\^ - A^) = Qg'{\^ - A^) = Rh'{\^ - A^) = const. (61) 

If /" = g" = h" = then the rank of the system also equals two. The requirement that the 
rank of the extended matrix equals two as well leads toc = 6 = a = 0. Thus, this case reduces 
to the previous one. 

Suppose now that f"g"h" ^ 0. Solving the linear system for P, Q, R we get 

_ 2{cf + bf + a)f _ 2(cff2 +bg + a)g' _ 2{ch^ + bh + a)h' 

(A2-A3)/" ' ^ (A3_Ai)5" ' (Ai-A2)/i" ■ 

Separating the variables in ()60p we ultimately obtain 

/'3 = ci52(/), g'^ = c2S\g), h'' = c3SHh), (62) 
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and 

(Ai - A^)(Ai - X^)S{f) (A2-Ai)(A^-A3)g(g) (A^-Ai)(A^-A^)5(^ 

2/' ' ^ 2g' ' 2/i' 

where is a polynomial of degree < 3, and q are arbitrary constants (the polynomial S{z) 
can be recovered from (A^ — A^)(A"^ — A^)(A^ — A^)S" = 6(cz^ + 6z + a)). Notice that the case 
(I6ip is a particular case of the above with S = const. 

We point out that the right hand sides of ()57p possess the following 5^(2, i?)-invariance, 



which can be used to bring the polynomial 5 to a canonical form. There are three cases to 
consider. 

Three distinct roots: in this case one can reduce 5 to a quadratic, S{x) = + g'3, so that 
the ODEs ([62]) imply / = p'(aiiti), g = p'{a2U2), h = p'iasus) where 27a? = 2ci and p is the 
Weierstrass p-function: (p')^ = ^p^ — 53. Up to a constant multiple, this leads to 

A* - A-^' p{aiUi)p{ajUj) _ ~ -^^ f ^ \ 



A - A- pyUiUiJp^UjUj) ^ A - A- L 

-"y ~ 77 A 77 V' ~ / V 2 — n^y^j'^j 



Qittj p'{aiUi) - p'{ajUj)' ^ \2 ^ ^ p'{aiUi) - p'{ajUj] 

The corresponding potential H{\\) is given by ()4ip . 

Double root: in this case one can assume S{x) = x, so that the ODEs (|62p imply / = 
(ai^i)'^, g = (021x2)'^, /i = (03^3)^, here 27af = Cj. Up to a constant multiple, this leads to 

^ _ (A'-A^Kn,- ^ _ (A^-A^>| 
The corresponding potential ^^^(u) is given by ()43p . 

Triple root: in this case S can be reduced to S = 1, so that the ODEs (j62p imply / = aiui, g = 
02^2, /i = a^us, here af = Ci. Up to a constant multiple, this leads to 

Hij = -. -, ifjj = - -, (63) 

aiaj[aiUi — ajUj) a-\aiUi — ajUj) 

and the corresponding potential H{u) is given by (j44p . Notice, however, that for this potential 
the expression (I52p equals zero. Formally, it should be considered as an example from the Case 
II below. 

Case II. This is the case when the expression (I52p equals zero, although both terms in ()52p are 
nonzero: 

i/l22i^233-/^113 = — ^/^112-f^223-f/^133 / 0; (64) 

an integrable example from this class is provided by 

\i _ \j 

Hi, = -, r; (65) 

aiaj[aiUi — ajUj) 

it appears in the 'triple root' case above. A detailed analysis below shows that this case possesses 
no other non-trivial solutions. Rewriting (I64p in the form 

ffl22 H233 Hiis _ _^ 
Hu2 H223 -f^l33 
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one can set 



Thus, 



H 



122 



H 



112 



H 



1 



233 



H' 



223 



H 



12 



1 



113 



133 



H 



1 
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where x = a{ui) — (3{u2), y = I3{u2) — 7(^3) and z = —x — y for some functions a, /3, 7 such that 
a' = l/l, (5' = 1/m, 7' = 1/n. Substituting ([66ll into (f32]) and integrating once, one gets 

(Ai - A2)p2 + (A3 _ Ai)i?2 



^11 

-f^22 
^^33 



(A2 - A3)PQi? 
(A3 - Ai)PQi? 

(A3 - \2)Q2 + (A3 _ Al)p2 

(Ai - X^)PQR 



+ Z/(mi,M3), 
+ r]{ui,U2). 



(67) 



Expressing six partial derivatives of the functions /i(u2,U3), z/(mi,M3), rj(ui,U2) from the six 
compatibility conditions djHa = diHij, and substituting them into the equations diJ = and 
(?2 J = 0, we obtain 



t(;ia'(ui) + W2-y'{u3) = 0, 
52(u^i)a'(ui) + d2iw2h'{u3) 



W3l3'{u2) + 1^47' (^^3) = 0, 
di{w-i)(3'{u2) + di{wi)i{u3) = ^, 



(68) 



where 



wi = (A^ - \^)Q{RP'Q' + QP'i?' - PQ'i?') 
it;3 = (A^ - X^)R{RP'Q' + QP'R' - PQ'R!) 



W2 = (A^ - \^)P{RP'Q' - QP'R' + PQ'R'), 
= (A^ - X^)P{RP'Q' - QP'R' - PQ'R'). 

The equations (f68]) 9 are obtained from ([68l) i upon differentiation by n2 and ui, respectively. 
Since a', (3' and 7' are nonzero, the system ()68p is consistent iff P, Q and R satisfy the following 
conditions: 

Wid2{w2) - W2d2{wi) = 0, wsdiiwi) - W4,di{w3) = 0. (69) 
Let us observe that the equations ()51|) (which also hold in this case) can be rewritten as follows: 



P_ 

P 



pq 

r 



rq r 
— + - 
p r 



pq 
r 



rp r 
— + - 
q r 



P_ 

p 



rp rq ^ q 



p 



(70) 



here we use the notation p = P'/P, q = Q' /Q, r = R'/R, and prime denotes the derivative 
of functions with respect to their arguments. Note that only two of the above equations are 
independent. Differentiating, for instance, the first two equations in (j70p by x and y and 
eliminating r and r' , one ends up at the following relations involving p and q: 



+ 



/\ 2 



2 {p' - q') 



(71) 
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Similarly, one can eliminate p and p' obtaining the analog of (|71|) for q and r. All these relations 
imply 

'■ij)-''-''iS)-'''-''iS)-''- '''' 

Thus, p, q and r must satisfy an ODE of the form 

//"-(/')'-/' + A; = 
where / = /(C) s-^id /c is an arbitrary constant. If /c = then 

/ = ^ 7 or / = - 

sm vQ C 

where is an arbitrary constant. Note that the second solution is a limit of the first as — > 0. 
If 7^ we have 

tanh/ciC ^ \ 

where sn is the Jacobi elliptic sine function: {sn')"^ = (l—sn?) (l — l/{ij^k)sn'^^. Using p = P'/P 
and integrating once, we obtain 

tanz/x sinh/cJj; n r cn 

P = C3X, P = C3 , -P = C3 J , or P = csdn - C3 . (73) 

Analogously, Q and can be obtained from these formulae by cycling the indices C3 — > ci ^ C2 
and the variables x — > y ^ z. Here cn and are the Jacobi elliptic functions cn{u^/kx; l/{u^k)) 
and dn{i>^/kx; l/(z^^/c)), and ci, C2 and C3 are arbitrary constants. It turns out that only linear 
and trigonometric solutions in (j73|) satisfy the condition (|69p . Thus, hyperbolic and elliptic 
solutions can be dropped. The substitution of the linear solution into one of the equations (f68l) i 
implies that the functions a, /3 and 7 must be linear. One recovers the solution (pH]) by setting 
ci = a2a3/(A^ — A^), C2 = aia3/(A^ — A^), C3 = aia2/(A^ — A^). Finally, the substitution 
of the trigonometric solution (|73p also implies that a, /5 and 7 must be linear, however, the 
compatibility conditions for the systems (f66]) and (fHTl) are not satisfied. 
Case III. This is the case when both terms in (|52p equal zero separately: 



9m2 c^?^3 5ui dui du2 dus 

Up to permutations of indices, we have to consider the following three subcases. 
Subcase 1: H12 = const 7^ 0. It follows from ([32]) that 

> 1 \ rr /x3 x2\Tr 

13- 



(A^ — ^^)Hi3H233 — {X^ — A^)f/'i2-f^223, (A^ — A^)i/23-f^l33 — (-^^ " A^)i/i2-f?l 



Differentiating the first equation with respect to ui we obtain -ffii3-f^233 = 0. If -?/233 = 
then the first equation implies -ff23 = const. Otherwise, it follows from the second equation 
that Hi3 = const. Without any loss of generality we assume that ^^23 = const ^ 0. Setting 
H12 = g(A^ — A^), 7^23 = p(A^ — A^), p,q = const, and substituting into ([5^ one arrives, up to 
the equivalence transformations, at the following potential H: 

H{ui,U2,U3) = {pui + qus) In {pui + qus) - ;^p(A^ - A^)(A^ - A^)uf- 

D 
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Aq(^X^ - X'){X^ - X^)ul +p{X^ - X^)U2U3 + q{X^ - X')uiU2. 
Subcase 2: H12 = f{ui), H23 = giu^). One can prove that in this case 

H{ui,U2,U3) = au\u2 + I3u2u\ + 7uf + 5u\ + u-^G 



for some constants a,/3,7, 5. The function G has to satisfy an equation of the form 

G"{x) = 

Z2 + ZsX-^ 

where Zi are some constants. If 23 = we have G{x) = . In this case 

a = (A2-A^), /3=(a2-A3), 7 = 0, <5 = ^(A2-A3)(A3-A1), 



which gives (I49p . The case 2:2 = is equivalent to the above. Otherwise, 

G{x) = {px + q) log (px + (?) + e(px + eg) log (px + eg) + e^(px + e^g) log (px + e^^ 
In this case 



« = (A2-Ai), fi = {X^-X\ J = r^^iX' - X'){X' - X% 6=r^{X^-X'){X'-X'), 



which gives (ISOj) . 

Subcase 3: i7i2 = f{ui), H13 = g{ui). A direct calculation shows that this case gives no 

non-trivial examples. 

This finishes the proof of Theorem 3. 

5.1 Dispersionless Lax pairs 

In this section we prove that the diagonalizability conditions (j32p imply the existence of the 
dispersionless Lax pairs (Theorem 4), and explicitly calculate Lax pairs for some of the most 
'symmetric' examples appearing in the classification list of Theorem 3. 

Example 1. Let us consider the quartic potential (j47p . 

H = {X^ - X^)ulul + (A^ - X^)ulul + (X^ - X^)ulu'l, 

which is a three-component generalization of the potential (jl9p from Theorem 1 (we have verified 
that this example possesses no natural four-component extensions). The corresponding system 
(j3ip has a Lax pair 



ipT = Xiai{C)ul + A2a2(0'"2 + ha3i^,)ul, = -ai{C)ul - a2{i)ul - «3(0^3; 

here ^ = ipY and the functions ai(^) satisfy the ODEs 

, 4ai , 4a2 , 4o3 

tti = 1-2, a2 = 1-2, 03 = h 2, aia2 + 0203 -|- 0301 = 0. 

03 «i 02 

Equivalently, 

4ai 4ai „ 2 10 

03 = — 0.2 = — r~r^' 2aiai = 3ai - 12. 

— 2 ai + 2 
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Without any loss of generality one can set 

ai = p(0> a2 = p(C + c), a3 = p(^-c) 

where p is the Weierstrass p-function, (p')^ = 4p'^ + 4, and c is the zero of p such that p(c) = 
0, p'(c) = 2. 

Example 2. Let us consider the potential (|44|) . 

H = — — ^ ^ (oitt^ — ajUj) ln(ainj — cijUj), 

which is a three-component generalization of the potential p8|) from Theorem 1. The corre- 
sponding system ([3TI1 possesses the Lax pair 



A* 1 

■0r = - ^ — ln(aini - tpy), V'x = ^ — ln(aini - ^y). 



This Lax pair appeared previously in [21j. 
Example 3. Let us consider the potential ()4ip . 



= - a^n ■). 

One can show that the corresponding system (I3ip has the Lax pair 

= - ^ ^f{aiUi,ijY), V^x = ^ \f{aiUi,il:Y) 
where the dependence of f{u,£,) on its arguments (here ^ = ^y) is governed by 

p'(n)-p'(0' P'(0-P'M 

Explicitly, one has 

1 e 
= - lna(u-0 + 7:lna(eu-0 + 7rlna(e\-0, 

DO b 

here a is the Weierstrass sigma- function: a' /o = (. In a different parametrization, this Lax pair 
appeared in [21] in the classification of dispersionless Lax pairs with movable singularities. We 
point out that both examples 2 and 3 generalize to n-component case in a straightforward way 
(allowing the summation to go from 1 to n). 

In fact, the following general result holds: 

Theorem 4 Any system (i3T]) satisfying the diagonalizability conditions [3^1 possesses a disper- 
sionless Lax pair. 
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Proof: 

We look for a Lax pair in the form (jlip . The compatibihty condition iptx = i^xt results in the 
following set of relations: 

/i = AV, /2 =X^92, f3 = X^93, (75) 

and 

fp9i = Hiigi + H2ig2 + -H'aiS's + 9pfi, 

fp92 = Hi2gi + H2292 + Hs293 + 9pf2, (76) 
fp93 = Hisgi + H2392 + Hssgs + gpf^, 

where we have set p = il^y, fi = dif, and gi = dig. The relations ([7^ and ([7^1) are equivalent 
to ()12p . Eliminating fp and gp from (I76p . one obtains a single algebraic constraint among the 
components gi,g2,93, which coincides with the left characteristic cone (j34p . The expressions for 
fp and gp obtained from the first two equations ()76p take the form 

„ _ (Hligi + gl2g2 + -^13^3) A^g2 - (Hugi + -g22ff2 + j^23g3) A^gl 

<7i52(A2-Ai) ' ^^^^ 

_ (-gjlgl + gl2g2 + -^13^3) g2 - (^^^12^1 + H2292 + i^23ff3) 9l 
~ 5152(A2 - Ai) 

Using the compatibility conditions fij = fji and /jp = fpi, we can express all second order 
derivatives of g in the form 

512 =513 = 523 = 0, 



511 



51 (^11151 + -H^ii252 + -^^11353) 
Hi2g2 + Hi^gs 



_ 92 {H22191 + H22292 + H22393) ^^g-^ 

Hl2gi + -^2353 



533 



51 (^^12351 + -f^22352 + -f^33253) (A^ - A^) 



+ 



^1352 (A3 - A2) + H2391 (Ai - A3) 

52 (-^^11351 + ^^12352 + ^^33153) (A^ - A3) 



^1352 (A3 - A2) + H23gi (Ai - A3) • 

It was already mentioned that the condition J = implies the decomposition of the left charac- 
teristic cone (pH) into a linear and quadratic factors, see (f35|) . We will assume that 51,52,53 lie 
on the quadratic branch, 

r = HisH23{X' - X^)9ig2 + Hi2H23{X^ - Ai)<7i53 + i^i2i^i3(A' - A3)5253 = 0. (79) 
One can verify that the differential consequences 

dT dT dT dT 

7^ = 0, 7^ = 0, 7^ = 0, Tr = 80 
oui 0U2 0U3 op 

hold identically modulo ([78|) , ([79]) and ([32]) . Finally, using computer algebra, it is straightforward 
to verify that the consistency conditions for the system ([TSll are satisfied identically modulo ([79ll 
and (|32p . This completes the proof of theorem [H 
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5.2 Hydrodynamic reductions 



The aim of this section is to prove that all examples listed in Theorem 3 possess infinitely 
many n-component hydrodynamic reductions parametrized by n arbitrary functions of a single 
variable. To do so one has to demonstrate the consistency of the relations dH), (jlOp where the 
characteristic speeds z/* and /i* satisfy the dispersion relation det{vl^ + + B) = Q, and SjU is 
the right eigenvector of the matrix i^^I^ + ^^A + B — see Sect. 2. 

Theorem 5 The diagonalizability conditions i32\) are necessary and sufficient for the existence 
of an infinity of n-component hydrodynamic reductions parametrized by n arbitrary functions of 
a single variable. 



Proof: 



The necessity follows from the general result of |T2] which states that, for a quasilinear system 
dH) , the diagonalizability is a necessary condition for the existence of an infinity of hydrodynamic 
reductions. 

The first step to demonstrate the sufficiency is to explicitly parametrize the dispersion curve 
Q, which we know to be a rational curve of degree three (see the Remark after Theorem 2). 
This can be done as follows. Let us first calculate the singular point UQ,fj,o on the dispersion 
curve. It corresponds to the situation when the rank of the matrix 1^/3 + fj,A + B drops to one. 
The associated left eigenvectors constitute a two-dimensional plane given by the first factor in 
the equation of the left characteristic cone (|35p . A simple calculation shows that i^q and fiQ can 
be obtained from the linear system 

1 H12H13 

1/0 + X flO = —rz Hii, 

2 H12H23 

UQ + X flQ = — H22, 

, ,3 H13H23 

fo + A /io = —[^ -H33; 

notice that these three relations are linearly dependent, indeed, multiplying the first by — A'^, 
the second by A^ — A^, the third by A^ — A^ and adding them together, one obtains J = 0, see 
()32p . Next, we parametrize the quadratic branch of the left characteristic cone ()35p in the form 

1 1 1 

91 = 7T1-. 77^' 92 = TT^- TT^, 93 



(A1+ 5)^23' (A2 + 5)ifi3' (A3 + 5)^12' 

here s is a parameter. The corresponding relation (|33p is equivalent to 

^ , X^+sHi2H23 , Al+ 5^13^23 . 

i^ + fiX +Hu + + -3- = 0, 

X^ + s X'^ + s H12 

^ X^^U , X-' + sHi2Hi3 , A^ + 5^13^23 . 

+ /xA +H22 + FT + T3- 7} = 0' 

X^ + s H23 A-^ + s H12 

^ \^^U , X^ + sHuHi3 , X' + sHi2H23 „ 

+ /xA^ + H33 + — — + = 0; 

X^ + s H23 X^ + s -H13 
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we point out that these three relations are also linearly dependent. Solving them for u{s) and 
fi{s) one obtains a rational parametrization of the dispersion curve: 





s 


H12H13 


s 


H12H23 


s 


H13H23 




+ s 


H23 


A2 + 


S Hi3 


A3 + 


s H12 




1 


H12H13 


1 


H12H23 


1 


H13H23 




+ s 


H23 


A2 + 


S Hi3 


A3 + 


s H12 



here fo and /^o £ire coordinates of the singular point. Thus, the characteristic speeds i^*(R) and 
/i*(R) can be represented in the form 

u\n) = u{s'), ^^\n) = ^,{s') (si) 

where s*, which are the parameter values of n points on the dispersion curve, are certain func- 
tions of the Riemann invariants: = s*(R). Since in our case the matrix 1^/3 + + B \s 
symmetric, the left characteristic cone coincides with the right characteristic cone. Thus, the 
right eigenvector corresponding to the point on the dispersion curve is 

1 1 1 ' 



(Ai + sOi/23' (A2 + s')/7i3' {\^ + si)Hi2 
and the relations ([8]) take the form 

A^ + H23 . . A^ + H23 . 
A^ + s* Hi3 \^ + Hi2 

Substituting (jSip into the commutativity conditions (jlOp and using (j82p one obtains the relations 

djs' = {...)djUi (83) 

^ 7^ J) where dots denote certain rational expression in s^,s^ whose coefficients depend on the 
second and third order derivatives of the potential H. For example, in the case of the quartic 
potential (j47p these relations take the form 

i _ 3(Ai + s'){X^ + ^^(A^ + s^)iX^ + gJ) 
" (Ai - A2)(Ai - A3) (si - ni 

By virtue of ([82]) and ([32]) . the consistency conditions djdiU2 = didjU2 and djdiUs = didjU3 
imply one and the same relation 

didjUi = {...)diUidjUi, (84) 

i ^ j, where, again, dots denote a rational expression in s*, whose coefficients depend on the 
second and third order derivatives of H. In the case (|47p . we have 

OiOnUi = OiUiOjUl, 

Ul 

where 

, . _ Qa^fi"^ + ki{a^(3 + a/j^) + ^2(0^ + 4a/3 + /g^) + ksja + (5) + 

(Ai-A2)(Ai-A3)(a-/3)2 
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fci = 3(A2 + + 2A^), k2 = (A^)2 + 2A^A2 + 2AU^ + X^X^, 

= 3A^(A^A2 + A^A^ + 2A2A^), ki = 6{X^fX^X\ 

The relations ([83]) and (fSl]) constitute the so-cahed Gibbons-Tsarev-type equations which govern 
hydrodynamic reductions of the system (jSip . The last step is to verify their consistency, namely, 
dkdjS^ = djdkS^ and didjdkUi = didkdjUi (without any loss of generality one can set i = 1, j = 
2, s = 3). If these consistency conditions are satisfied identically, the system ()83p . (j84p will be in 
involution, with the general solution depending on 2n arbitrary functions of a single variable. Up 
to reparametrizations i?* P{R^) this gives an infinity of hydrodynamic reductions depending 
on n arbitrary functions. 

We have verified the consistency for all examples appearing in Theorem 3. In fact, rather 
than considering them case-by-case, one can give a unified proof of the consistency using only the 
diagonalizability conditions (f32|) . To do so one needs to bring the system (|32]) into a passive form. 
It turns out that all higher order partial derivatives of the potential H can be expressed in terms 
of the second order derivatives Hij and the 4 third order derivatives, say, Hi22-, -ffiis, -^223; -^233- 
Second order derivatives are constrained by a single algebraic equation J = 0, while the values 
of H and its first order derivatives Hi are arbitrary. This calculation shows that the generic 
solution of the system (j32p should depend on 13 arbitrary constants, which is in full accordance 
with the results of Section 5. The computation of the expressions ()83p and (I84p . as well as 
the verification of the consistency conditions have been performed modulo this passive form. 
This means that all partial derivatives of H except the basic ones were eliminated, and the 
basic derivatives were considered as independent variables related by a single algebraic equation 
J = 0. An intense computer calculation shows that all compatibility conditions are identities in 
the basic derivatives. 

6 Hamiltonian systems in 3 + 1 dimensions 

In this section we establish a number of non-existence results for integrable Hamiltonian systems 
of hydrodynamic type in 3-|-l dimensions. We will begin with a two-component case. According 
to the results of [19], there exists a unique two-component Hamiltonian operator of hydrodynamic 
type which is essentially three-dimensional. Up to a linear transformation of the independent 
variables it can be cast into a canonical form 

/ d/dx \ / d/dz \ 
\ d/dy d/dz J ■ 

The corresponding Hamiltonian systems -|- P{h^^) = take the form 

ul + + {h2)z = 0, u^t + {h2)y + (hi), = 0. 

Applying the Legendre transform, ui = hi, U2 = h2, H = v}hi + v?h2 — h, one can rewrite 
these equations in the equivalent form 

{ui)x + {U2)z + {Hi)t = 0, {U2)y + (ni), + {H2)t = 0. (85) 

Notice that H{ui,U2) is defined up to an arbitrary quadratic form (all quadratic terms in H 
can be eliminated by appropriate linear changes of the independent variables). Our first result 
is the following 
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Theorem 6 Any integrable system is necessarily linear (that is, the potential H is quadratic 
in ni, M2 j. 



Proof: 

Our strategy will be to consider reductions of the system (|85|) to various (2 + l)-dimensional 
systems. In fact, it will be sufficient to look at reductions governing traveling wave solutions. If 
the original system (|85p is integrable, all such reductions must be integrable as well. Since the 
integrability conditions for (2 + l)-dimensional two-component systems of hydrodynamic type 
are explicitly known [llj, this will provide a set of necessary conditions for the integrability 
of the system ([85]) . It turns out that these conditions are very strong indeed, leading to the 
non-existence of non-quadratic integrable potentials H. 

Setting in the equations ()85|) dz = fJ-dt, which is equivalent to seeking solutions in the form 
u(x,y,t + fj,z), one obtains a (2 -|- l)-dimensional Hamiltonian system 

{Ui)x + (Hi + fJ.U2)t = 0, {U2)y + {H2 + iiui)t = 0, 

with the Hamiltonian density H(ui,U2) + fJ-uiU2. According to our philosophy we have to require 
that it is integrable for an arbitrary value of the parameter fi. The integrability conditions ()15p 
readily imply that the corresponding H must be cubic in ui,U2, and Theorem 1 tells us that 
the only two 'suspicious' cases to consider are H = and H = ^niM^ (recall that we ignore 
quadratic terms in H). In the first case the system (j85p takes the form 

iui)x + {U2)z = 0, {U2)y + {Ui)z + U2{u2)t = 0. 

Setting here x = y (this amounts to seeking traveling wave solutions in the form u(x -|- y, z, y), 
one obtains a (2 -|- l)-dimensional system 

{Ui)x + {U2)z = 0, {U2)x + {Ul)z + U2{u2)t = 0. (86) 

We recall that the paper [llj provides a complete set of the integrability conditions for two- 
component hydrodynamic type systems represented in the form 




The integrability conditions constitute a complicated over-determined system of PDEs for the 
coefficients a, q, r, s as functions of v, w. Representing the equations (f86]) in the form 




one can verify that these integrability conditions are not satisfied. Thus, the (3-|- l)-dimensional 
system corresponding to H = is not integrable. Similarly, for H = ^uiu^ the system ([85]) 
takes the form 

{Ul)x + {U2)z + U2{u2)t = 0, {U2)y + {ui)z + U2{ui)t + Ui{u2)t = 0. 
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Setting, again, x = y, and changing to the new dependent variables v = ui + U2, w = U2 — ui, 
one obtains the system 



3v+w 




v—w 
4 

_ v+3w 
4 



)l 


f - \ 




\ w ] 



which also does not satisfy the integrability conditions. This finishes the proof of Theorem 4. 

Our next result shows that any three-component (3 + l)-dimensional integrable Hamiltonian 
system associated with a non-singular Poisson bracket of hydrodynamic type is either linear or 
reducible. Any such system can be brought to a canonical form 



u\ + (/ii), 
with the Hamiltonian operator 



0, ul + {h2)y = 0, + (/13), 



d/dx 
d/dy 
d/dz 



0, 



(87) 



Performing the Legendre transform one obtains 

{Hi)t + (ni)^. = 0, {H2)t + {U2)y 

or, in matrix form, 



0, 



where the 3x3 matrices Ai are given by 



= 0, 



0, 



An 





H12 


Hi3 


H12 


H22 


H23 


Hi3 


H23 


H33 






Theorem 7 Any integrable (3 -|- 1)- dimensional Hamiltonian system (81) is either linear or 
reducible. 



Proof: 

As a necessary condition for integrability, one has to require the vanishing of the Haantjes tensor 
for an arbitrary matrix of the form 

(^0 + A^i + f3A2 + 7^3)'^ (^0 + AAi + M2 + 7^3), 

which is equivalent to the vanishing of the Haantjes tensor for any matrix A(Ao -|- A) where A 
and A are arbitrary 3x3 constant coefficient diagonal matrices. Computing the Haantjes tensor 
end equating to zero coefficients at different monomials in the diagonal entries of A and A, one 
obtains that either all third order derivatives Hijk are identically zero (this corresponds to linear 
systems), or Hij = Hik = for some i ^ j ^ k (this corresponds to the reducible case). 

We would like to conclude this section by formulating the following general 
Conjecture There exists no non-trivial integrable Hamiltonian systems of hydrodynamic type 
in 3 -|- 1 dimensions corresponding to a local Poisson bracket of hydrodynamic type and a local 
Hamiltonian density. 
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7 Concluding remarks 



We have found a broad class of non-trivial potentials leading to integrable Hamiltonian systems 
of hydrodynamic type in 2 + 1 dimensions. There is a number of natural problems arising in 
this context, in particular: 

— Describe the structure of the corresponding Hamiltonian hierarchies. The main difficulty 
here is the non-locality of higher symmetries / conservation laws. 

— Construct the associated Hamiltonian hydrodynamic chains. This requires the introduction 
of a canonical set of non-local variables reducing all higher flows of the hierarchy to infinite- 
component systems of hydrodynamic type. 

— Construct dispersive deformations of the examples arising in the classification, especially 
those with 'elliptic' Lax pairs. 

— Study the behavior of exact solutions coming from hydrodynamic reductions. 
We hope to address some of these questions elsewhere. 
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